Scaling of Cluster and Backbone Mass Between Two Lines in 3d Percolation 
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We consider the cluster and backbone mass distributions between two lines of arbitrary orien- 
tations and lengths in porous media in three dimensions, and model the porous media by bond 
percolation at the percolation threshold Pc- We observe that for many geometrical configurations 
the mass probability distribution presents power law behavior. We determine how the characteristic 
mass of the distribution scales with such geometrical parameters as the line length, w, the minimal 
distance between lines, r, and the angle between the lines, 9. The fractal dimensions of both the 
cluster and backbone mass are independent of w, r, and 9. The slope of the power law regime of 
the cluster mass is unaffected by changes in these three variables. However, the slope of the power 
law regime of the backbone mass distribution is dependent upon 9. The characteristic mass of the 
cluster also depends upon 9, but the characteristic backbone mass is only weakly affected by 9. We 
propose new scaling functions that reproduce the 9 dependence of the characteristic mass found in 
the simulations. 



I. INTRODUCTION 



Since the 1950s, the percolation model has been ap- 
plied to many disordered systems |l|-^ , and continues to 
be useful today. Here we use percolation theory to ana- 
lyze the cluster and the backbone mass distributions of 
clusters that are connected in configurations of the type 
shown in Fig. l|, configurations in which the two lines 
are connected by occupied bonds. The cluster mass is 
the number of bonds that are connected to the two lines 
(Fig. |l|a). The backbone of the cluster is the set of bonds 
that are connected to the two lines through independent 
paths (i.e., paths that have no common bond [^||). For 
configurations of two points, the distributions of vari- 
ous quantities have been studied [p|-p^. Recently the 
distribution of the shortest paths between two lines has 
been studied for a three-dimensional cubic lattice [ p^ , 
and here we calculate the mass and backbone distribu- 
tions. 

The motivation for this study is its relevance to tech- 
niques of oil recovery in oil fields |16|. A common tech- 
nique used in oil recovery is the injection of fluid into 
the ground at one site in the field in order to force oil 
out of the ground at another site nearby (Fig. |l|). It is 
common to inject the fluid along a portion of the length 
of the injection well and to collect the oil along a por- 
tion of the length of the production well (as opposed to 
injecting and collecting at single points on the wells). In 
our model, each line represents a well in the oil field. 
One line represents the injection well, and the other the 
production well. In many cases the oil reservoir is ex- 
tremely heterogeneous, and the percolation model is ap- 
propriate. Separation of the rocks into two types — high 
permeability ( "good rock" ) and low or zero permeability 



("bad rock") — can be accomplished at the outset, with 
the good rock represented by occupied bonds and the bad 
rock represented by unoccupied bonds. The connected 
mass represents the total oil in the reservoir connected 
to the two wells, and the backbone mass the recoverable 
oil. 



II. SIMULATIONS 

We perform a numerical study of the system using 
Monte Carlo simulations. We specify two sets of points 
representing lines in a simple cubic lattice to be the wells 
and we grow the cluster from these two lines of seeds. If 
the growth of either cluster stops before the two clusters 
connect, we discard the realization. For realizations in 
which the two clusters connect, the simulation ends either 
when the cluster growth stops naturally, or when the clus- 
ter mass reaches some specified limit. To eliminate finite 
size effects, we use the techniques of Ref. Q to simulate 
systems on lattices of large enough size that the clus- 
ters never reach the edge of the lattice. We perform the 
simulations at the percolation threshold, Pc = 0.2488126 
|]l7| . The configurations are characterized by three pa- 
rameters: length w, angle 9, and minimal distance r [see 
Fig. |^(a)]. For each configuration, we run at least 10^ 
non-discarded realizations. We calculate the cluster and 
backbone mass for each of these realizations as exempli- 
fied in Fig. pl( 



III. GENERAL OBSERVATIONS 

In both the cluster mass and backbone mass distribu- 
tions we expect to observe an initial cutoff due to the 
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fact that these masses cannot be smaller than the dis- 
tance r. In the backbone distributions, we expect to 
see a second cutoff due to the fact that the backbone 
mass cannot be greater than the cluster mass at which 
we stop the simulations. In addition, for both types of 
distributions we expect to observe a regime that exhibits 
power-law behavior. These general features of the distri- 
butions have been observed in the distributions for other 
quantities [p[-p^. The quantities of interest are (i) the 
most-probable value of the distribution (the maximum), 
the scaling of which will be determined by the fractal di- 
mensions of the quantities measured, and (ii) the slope of 
the power-law regime. For clusters grown from a single 
point, the slopes of the power law regimes of the cluster 
and backbone distributions are r — 1 and tb — 1, where 
r is the Fisher exponent and tb the corresponding expo- 
nent for the backbone. The fractaljdimensions and power 
law regime slopes are related by [ 



T - 1 = 



TB-l 



d 
d_ 



(1) 
(2) 



where d is the dimension of the system, and df and ds are 
the fractal dimensions of the cluster and the backbone, 
respectively. For d = 2>, estimates for these exponents 
are |l8|jl|] 



a/ 
ds 

T- 1 
TB-l 



2.524 ±0.008 
1.855 ±0.015 
1.189 ±0.004 
1.617 ±0.013. 



(3a) 
(3b) 
(3c) 
(3d) 



IV. CLUSTER MASS 



A. Parallel Wells 



We perform the same analysis as before, and ob- 
tain similar results [see Fig. §(a)], i.e., power law 
distribution for P(m\w) with a slope « —1.18 and 
fractal dimension « 2.55 [Fig. §(b)]. 

We now study cases intermediate to those studied in 
(i) and (ii). In Fig. 5a, we plot the distribution of cluster 
mass for configurations in which r = 16 and we vary w 
from to 64. For small w, the distributions are essen- 
tially unchanged, but for w ^ r, the distributions scale 
with the exponent df (Fig. 5a and 5b). 

We now develop a scaling form for the dependence of 
the characteristic mass m* on r and L, the system size. 
Without loss of generality we can write 



(4) 



This form is consistent with the scaling of the cluster 
mass. That is, if 



r' = ar 
w' = aw 



(5) 



then 

m*[r',L') = 



/(-) 



= a'^fm*{r, L). 
(6) 

We further assume that m*{r, L) can be written as 



m*{r, L) = 
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(7) 



since w and r become irrelevant variables for r ^ w and 
r ^ w, respectively. Thus we expect 
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r ^ w 

constant r w 



(8) 



In order to gain insight into the general behavior, we 
first study parallel wells (6* — 0) see Fig. ^(a). We con- 
sider first the following limiting cases: 

(i) w r — In this case we approximate the configu- 
ration by two points (see Fig. ^(b)). In Fig. |^(a) we 
show the mass probability distribution P{m\r) for 
w; = and r = 1, 2, 4, 8, 16, 32 and 64. The distri- 
bution shows a maximum followed by a power-law 
regime with slope —1.18, consistent with Eq. (|3^). 
We study also how the characteristic mass m* , cor- 
responding to peak of the distribution, scales with 
the distance r. The log-log plot of m* vs r in Fig. 3b 
indicates that m* scales with exponent ds = 2.6 
which is consistent with Eq. (|3a]). 

(ii) w ^ r — For this case [see Fig. ||(c)] we approx- 
imate the configurations by r = (a single line). 



B. Non-Parallel Wells 

We now study non-parallel wells. The results for 
the mass probability distribution P(rn\0) are shown in 
Fig. ^(a). We find that the power law regime is consis- 
tent with a slope —1.18 independent of 6. We consider 
also the peak of these distributions, analyzing how m* 
evolves with 0. Figure ^(b) shows the dependence of m* 
versus 6. In this case we do not find a power distribu- 
tion for the scaling behavior m* versus 9. m* increases 
rapidly for small values of 6, and for larger asymptoti- 
cally approaches a limiting value at = vr. 

We now suggest a functional form for the dependence 
of the characteristic mass m* on 9. Without loss of gen- 
erality we can write 

m*{9)=m*i0)[f{9)]''f. (9) 
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Since the configuration for = tt is simply a single 
straight line twice the length of the single line for 6* = 0, 
we expect 



(10) 



(11) 



TO*(7r) = m*{0)2'^f. 
We then are motivated to write 

m*{9)=m*{0)[l + g{e)]'^f, 
where g{9) is monotonic and 

5(0) = 

and 

9{^) = 1- 



A first guess at a functional form for g(6) is some power 
of sin(6'/2) but no power seems to fit the data of Fig. ||(b) 
well. A functional form that fits better is 



h{0) 



^ sm , „ 
2 V 2 



The final functional form for m*{9) is thus 
m*{e) = m*(0) I 1 + sin 



— sm , 
2 V 2 



0.4N 



(12) 



(13) 



where the exponent 0.4 is obtained by the power law fit in 
Fig. |^(a). We note that there is no a priori justification 
for this form; it simply satisfies the appropriate bound- 
ary conditions and fits the simulation results reasonably 
well, as shown in Fig. 0(b) . 



V. BACKBONE MASS 

We proceed to analyze the scaling of the backbone 
mass distributions as we did for the cluster mass. 



A. Parallel Wells 

(i) w <^ r — We present in Fig. P(a) the back- 
bone probability distribution P{mB\r). As seen in 
Fig. ^(a) and (b), the power law regimes of these 
plots are consistent with Eq. (3d), and the distri- 
butions scale with an exponent consistent with cIb- 
Similar results have been found previously jl^ ] but 
are shown here for completeness. 



(ii) w ^ r — We show in Fig. ||(a) the backbone prob- 
ability distribution P{mB\'w) for this case. Again, 
the power law regimes of these plots are consistent 
with Eq. ( ^ ) and the distributions scale asymp- 
totically with an exponent consistent with ds 
[Fig. ||(b)] although there are large corrections-to- 
scaling as seen in the small w behavior of Fig. ^(b). 



Figures ^(a) and ^(b) are also indicative of the behav- 
ior of cases intermediate to the limiting cases of (i) and 
(ii) above. For w <^ r, changing w will have essentially no 
effect on the distributions; for w ^ r, the distributions 
will scale with ds as w is changed. 



B. Non-Parallel Wells 



We show the results for the backbone probability dis- 
tribution P{mB\0) in Fig. |l^(a) for various values of 9. 
In contrast to the distributions of cluster mass, the dis- 
tributions of backbone mass exhibit power-law regimes, 
the exponents of which depend on 9 while the character- 
istic mass m*, is essentially constant as a function of 9. 
Thus an appropriate functional form for P{mB\9) is 



(14) 



where /i and /2 are cutoff functions. The first cutoff 
function, /i, reflects the fact that the backbone mass 
must always be at least equal to the distance r between 
the two points; the second cutoff function, /2, reffects the 
fact that the backbone mass is bounded because of the 
finite size, L, of the system. Similar behavior has been 
observed for the distributions of shortest paths between 
two fines Q. 

In order to determine the varying slope more accu- 
rately we perform simulations for various 9 for r = 1, 
which results in the largest power-law regime. The re- 
sults of these simulations are shown in Fig. |o|(b). In 
Fig. |o|(c) we plot the power-law regime exponent, gs 
vs 9. For 9 = 0, the configuration is that of parallel 
lines, and the exponent is tb- For 9 = tt, the exponent 
decreases to a value of about 0.84. The marked differ- 
ence in these two exponents is shown clearly in Fig. 11, 
in which we plot for fixed r = 1, 9 = 180°, P{mB\r) for 
various values of w. The larger the value of w, the later a 
crossover occurs from behavior reflecting a configuration 
of 2 lines with 9 = 180° to a configuration effectively of 
2 points, with power-law-regime exponent — 1. 



VI. DISCUSSION 



We have analyzed the distributions of cluster mass and 
backbone mass for various configurations of 2-line 3d per- 
colation clusters. The behavior of the cluster mass distri- 
butions is not remarkable. On the other hand, we have 
found that the exponent of the power-law-regime for the 
backbone mass distributions is dependent on the angle 9 
between the lines. It remains to develop a theory which 
can predict the specific dependence of this exponent on 
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(b) 



FIG. 1. (a) Illustration of well geometry, (b) Examples of a percolation cluster with two line wells with parameters r = 2, 
6 = 90°, and w = \/50. The filled sites are members of the percolation cluster, which has a mass of 52. Solid lines form the 
backbone, which has a mass of 28. 
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FIG. 2. Parallel well examples, (a) general case (b) w <C r (c) w S> r. 
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FIG. 3. (a) Cluster mass distribution P{m\r), for the percolation cluster for the two points case {w = 0), for several values 
of r (r = 1, 2, 4, 8, 16, 32, 64). The line of slope —1.18 denotes the theoretical expectation, (b) Scaling behavior of m*, the most 
probable mass, as a function of r. The fitted slope of 2.61 is consistent with the fractal dimension df = 2.54 
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FIG. 4. (a) Cluster mass distribution P{m\w) for (r = 0) for several values of w (w = 1, 2, 4, 8, 16, 32, 64, 128, 256, 512). (b) 
Scaling behavior of m* as a function of w. 
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FIG. 6. (a) Mass distribution P{m\0) for r — and w — 32 for the percolation cluster to the general non parallel wells for 
several values of the angle 9. (b) Corresponding scaling behavior of m* versus 9. 
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FIG. 7. (a) Determination of exponent in Eq. (12), where x = sin[(7r/2) sm{9/2)]. (b) Comparison of functional form rn{6) 
(solid line) versus observed data. 
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FIG. 8. (a) Backbone mass distribution P(ms|r) for the percolation cluster for two points and r = 1,2,4,8,16,32,64,128 
(b) Scaling behavior of mg versus r. 
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FIG. 9. (a) Backbone mass distribution P{mB\w) for parallel wells for r = 1 and w = 1,2,4,8,16,32,64,128. (b) The 
corresponding scaling behavior of rrig versus w. 
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FIG. 10. (a) Backbone mass distribution P(ms|S) for non-parallel wells with r = 8, w = 64, and several values of 6. The 
cutoff of cluster growth is at a cluster mass of 2^**. (b) Backbone mass distribution P{mB\0) for the case of nonparallel wells 
with r = 1, w = 64 for various 9. The cutoff of cluster growth is at a cluster mass of 2^". (c) Power law exponent gB{0), defined 
in Eq. (H), for the corresponding backbone mass distribution presented in (b). 
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FIG. 11. Backbone mass distribution P(ms|w) for the non parallel wells for fixed 6 (= 180°) and several values of w. 
The larger the value of w, the later a crossover occurs from behavior reflecting a configuration of 2 lines with 6 = 180° to a 
configuration effectively of 2 points, with power- law-regime exponent rs — 1. 
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